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Abstract. For a polygon P with n vertices, the vertex guarding problem 
asks for the minimum subset G of P's vertices such that every point in 
P is seen by at least one point in G. This problem is NP-complete and 
APX-hard. The first approximation algorithm (Ghosh, 1987) involves 
decomposing P into cells that are equivalence classes for visibility 

from the vertices of P. This discretized problem can then be treated as 
an instance of set cover and solved in C(w^) time with a greedy 0(log n)- 
approximation algorithm. Ghosh (2010) recently revisited the algorithm, 
noting that minimum visibility decompositions for simple polygons (Bose 
et at, 2000) have only 0(n'^^ cells, improving the running time of the 
algorithm to for simple polygons. 

In this paper we show that, since minimum visibility decompositions 
for simple polygons have only C(w^) cells of minimal visibility (Bose et 
al., 2000), the running time of the algorithm can be further improved 
to O(n^). This result was obtained independently by Jang and Kwon 
(2011). We extend the result of Bose et al. to polygons with holes, show- 
ing that a minimum visibility decomposition of a polygon with h holes 
has only 0(^{h + l)j'i'^) cells and only 0((/i + l)'^"'^) cells of minimal 
visibility. We exploit this result to obtain a faster algorithm for vertex 
guarding polygons with holes. We then show that, in the same time com- 
plexity, we can attain approximation factors of ©(log log opt) for simple 
polygons and C((l + log {h + 1)) log opt) for polygons with holes. 



1 



1 Introduction 

Art gallery problems, i.e., polygon guarding problems, are motivated by the 
question, "How many security cameras are required to guard an art gallery?" 
The art gallery is modeled as a connected polygon P. A camera, which we 
henceforth call a guard, is modeled as a point in the polygon, and we say that a 
guard g sees a point q in the polygon if the line segment Ijq is contained in P. 
The visibility polygon of a point p, denoted Vis(p), is the set of points in P that 
see p. We call a set G of points a guarding set if every point in P is seen by some 
g G G, i.e., if UgeG ^is(5) — ^- Let V{P) denote the vertex set of P and let dP 
denote the boundary of P. We assume that P is closed and non-degenerate so 
that V{P) CdPcP. 

We consider the minimization problem that asks, given an input polygon 
P with n vertices, for a minimum guarding set for P. Variants of this problem 
typically differ based on what points in P must be guarded and where guards can 
be placed, as well as whether P is simple or contains holes. Typically we want to 
guard either P or dP, and our set of potential guards is typically V{P) (vertex 
guards), dP (perimeter guards), or P (point guards). This paper concerns the 
variant in which we must guard all of P from vertices of P. For results on art 
gallery problems not related to minimization problems we direct the reader to 
O'Rourke's book [22], which is available for free online. 

1.1 Related Work 

Hardness Results The problem was proved to be NP-complete for polygons 
with holes by O'Rourke and Supowit ^23]. For guarding simple polygons it was 
proved to be NP-complete for vertex guards by Lee and Lin ^21, ; their proof was 
generalized to work for point guards by Aggarwal [1 . This raises the question of 
approximability. There are two major hardness results. First, for guarding sim- 
ple polygons, Eidenbenz [TU] proved that the problem is APX-complete, mean- 
ing that we cannot do better than a constant-factor approximation algorithm 
in polynomial time unless P = NP. Subsequently, for guarding polygons with 
holes, Eidenbenz et al. proved that the minimization problem is as hard to 
approximate as set cover in general if there is no restriction on the number of 
holes. It therefore follows from results about the inapproximability of set cover 
[2511312) that, for polygons with an unbounded number of holes, it is NP-hard 
to find a guarding set of size o(log n) . These hardness results hold whether we 
are dealing with vertex guards, perimeter guards, or point guards. 

Approximation Algorithms Ghosh [14 provided an 0(log 7i)-approximation 
algorithm for guarding polygons with or without holes with vertex guards. His 
algorithm decomposes the input polygon into a polynomial number of cells such 
that each point in a given cell is seen by the same set of vertices. This discretiza- 
tion allows the guarding problem to be treated as an instance of set cover and 
solved using general techniques. In fact, applying methods for set cover devel- 
oped after Ghosh's algorithm, it is easy to obtain an approximation factor of 
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O(logOPT) for vertex guarding simple polygons or 0(log ft, log opt) for vertex 
guarding a polygon with h holes. 

When considering point guards or perimeter guards, discretization is far more 
complicated since two distinct points will not typically be seen by the same set 
of potential guards even if they are very close to each other. Deshpande et al. 
[S] obtain an approximation factor of O (log opt) for point guards or perimeter 
guards by developing a sophisticated discretization method that runs in pseu- 
dopolynomial timcQ Efrat and Har-Peled [3] provided a randomized algorithm 
with the same approximation ratio that runs in fully polynomial expected time; 
their discretization technique involves only considering guards that lie on the 
points of a very fine grid. 

1.2 Range Spaces and Discretization 

Guarding problems can naturally be expressed as instances of set cover or hitting 
set. We wish to model an instance of a guarding problem as an instance of hitting 
set on a range space S = {X, TZ), constructed as follows. X is equal to the set Sg 
of potential guard locations. For each point p that needs to be guarded, Rp is the 
set of potential guards that see p. Now TZ = {Rp : p g St}, where St is the set 
of points that must be guarded. For the vertex guarding problem, Sg = V{P) 
and St = P. 

We assume Sg is finite; in our case \Sg\ — n. If St is not finite, e.g., when 
St — P, we need to discretize it. The goal of discretization is to find a finite 
representative subset S'j, C St such that any subset of Sg that guards 5'^ also 
guards St- With such a set we are able to forget about St and focus on the 
finite range space {Sg, {Rp ■ P G S^}) induced by Sg and S^. 

We consider general techniques for solving hitting set for finite range spaces. 
For a finite range space S = {X, TZ) the time complexity typically depends on 
\X\ (i.e., the number of elements) and \TZ\ {i.e., the number of ranges). 

1.3 Visibility Decompositions for Polygons 

For points p,q ^ St, we say that p and q are equivalent if and only if Rp = Rq. 
In a decomposition of a polygon into cells, we say that a cell is an equivalence 
cell if all points are equivalent. A natural discretization strategy is to partition 
St into a finite number of sets that are closed under equivalence, and then to 
build a subset S'^p by taking one representative point from each set. A subset of 
Sg guards St if and only if it guards S'rp. 

Ghosh [14] did this for the vertex guarding problem in which Sg — V and 
St = P. His algorithm decomposes the input polygon into a polynomial number 
of cells such that each point in a given cell is seen by the same set of vertices. For 
two distinct vertices that see each other, consider the line through them. The 
set of all such lines decomposes P into a number of equivalence cells (see Figure 



^ It is a pseudopolynomial-time algorithm in that its running time may be linear in 
the ratio between the longest and shortest distances between two vertices. 
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[I]). Using a simple inequality for general line arrangements it can be seen that 
the number of cells is 0(n*). Ghosh originally used his discretization technique, 
along with the greedy set cover approximation algorithm (see, e.g., '29' pp. 16- 
19]), to provide a ©(log n)-approximation algorithm. This decomposition can be 
generalized to work for any finite set Sq of potential guards — simply shoot a 
ray from every point p e So through every vertex seen by p. In each of the 
0(|5Gp?^^) resulting cells, any two points see the same subset of Sg- 



Fig. 1. A polygon decomposed into 
O^n*) cells by cutting along any line 
passing through two vertices that see 
each other. The shaded region is a hole. 



Fig. 2. A polygon decomposed into the 
0[n^) cells of its minimum visibility de- 
composition T>v{P,V{P)). The shaded 
region is a hole. 



Bose et al. IS] introduced a minimum decomposition with fewer cells (some of 
their results were obtained independently by Guibas et al. [TB], but we focus on 
the results as stated and proved by Bose et al.). Let the visibility decomposition, 
denoted 'Dv{P, So), be the minimum decomposition of a polygon into equiva- 
lence cells with regard to Sq. It is minimum in that the union of cell boundaries 
in this decomposition is exactly equal to UpeSc ^(^^^(p))- ^he decomposition 
can be constructed as follows. For a point p G Sq that sees a vertex v, instead 
of cutting along the entire ray shot from p through v, we leave the line segment 
pv and cut only from v until we hit dP (see Figure [2]). Bose et al. call such a cut 
a window of point p and they note that the boundary of Vis(p) consists only of 
windows of p and parts of dP. They proved an upper bound of 0{n^) for the 
number of cells in T>y{P, V{P)) for a simple polygon P. 

Ghosh |15j recently revisited his algorithm, using these minimum visibility 
decompositions to improve the running time. His updated algorithm guarantees 
an approximation factor of ©(logn) and runs in 0(ji'^) for simple polygons and 
0(n^) for polygons with holes. 
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1.4 Novel Contributions 

In this paper we exploit another result of Bose et al. [51, namely the fact that, 
while minimum visibility decompositions for simple polygons can have 
cells, they have only ©(n^) cells of minimal visibility. We use this to further 
improve the running time of Ghosh's algorithm from ©(n^) to 0(n'^). 

We then extend the result of Bose et al. to polygons with h holes, param- 
eterizing bounds not only by n but also by h. We show that a minimum visi- 
bility decomposition of a polygon with /i > holes has only 0{{h -I- l)n?') cells 
(Theorem[l]) and only 0{{h + 1)^ n^) cells of minimal visibility (Theorem 2L We 
exploit this result to improve the running time of Ghosh's algorithm from ©(n^) 
to 0{{h -\- l)^n^); this is a strict asymptotic improvement for h — o[n). 

Having presented algorithms with faster running times, we turn our attention 
to improving the approximation factor. We show that, with the same time com- 
plexity bound of ©((/i -I- l)^n^), we can apply standard random sampling tech- 
niques for range spaces of bounded VC dimension to obtain an approximation 
factor of ©((1 -I- log (/i -t- 1)) log opt) for polygons with h > Q holes (Theorem 

HI- 

Finally, we show that, with the same time bound of 0(n'^), we can achieve 
an approximation ratio of ©(log log opt) for simple polygons using the improved 
£-net finders of King and Kirkpatrick [T3] (Theorem [6]) . 

The author originally suggested the idea of exploiting the number of sinks to 
improve running time in a recent thesis [18] . This idea was also used indepen- 
dently by Jang and Kwon pTj. who obtained the same time complexity as us 
for simple polygons. Jang and Kwon also consider the problem of edge guards, 
whereas we do not. However, they do not consider improvement of the ©(logn) 
approximation factor or consider polygons with holes. 

1.5 Model of Computation 

We assume the real- RAM model of computation [21] . We also assume that the 
polygon and any holes are non-degenerate, and that vertices of the polygon are 
in general position, i.e., no three are coUinear. 

2 Polygon Decompositions 
2.1 Notation 

For a polygon P with h > holes, P U dP is the complement of the relative 
interior of P, or the closure of the complement of P. We use Cq, . . . , to denote 
the h + 1 components of P U dP. Cq is the polygon's exterior and Ci, . . . ,Ch are 
the holes of the polygon. 

In addition to left windows, right windows, left pockets, and right pockets 
defined by Bose et al. for simple polygons [5] (see Figure [s]), we have T- windows. 
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Fig. 3. Windows and pockets of a vertex v. Windows are represented by coloured line 
segments: left windows are blue, right windows are green, and T-windows are red. Left 
and right pockets are shaded blue and green respectively. 

which are trans-component windows. A window is a T-window if its two end- 
points are on different components. There are no pockets associated with T- 
windows. 

A window of a vertex v has two endpoints; the endpoint closer to v is called 
the base and the endpoint farther from v is called the end. The window half 
plane of a left (resp. right) window of v whose base is h is the half plane on the 
left (resp. right) side of the oriented line from v to b. 

2.2 Bounding the Number of Cells 

Here we extend the results of Bose et al. [5] concerning the number of cells in 
a visibility decomposition. Our proofs are, to a great degree, extensions of their 
proofs. 

Lemma 1 ( [Sj ) • No point z in a pocket Q of a window w is visible to any point 
y that is inside w's half plane but outside Q. 

Proof. The proof of (Sj Lemma 1] works without modification when generalizing 
to polygons with holes. For completeness we restate it here. 

Let w have b as its base and e as its end. Since QU w is a polygon, the line 
segment 'zy intersects QUw at least once. As yz cannot intersect be, it intersects 
some other line segment on the boundary of the polygon. So, y and z are not 
visible with respect to the chain Q. 

Lemma 2. For any vertex v there are at most 2h T-windows of v. 
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Proof. Consider the component sequence of a vertex v obtained as follows. Let 
r be a ray emanating from v that bisects the external angle at v (or, if w is a 
vertex of a hole, the angle interior to the hole). From its starting position, the 
r rotates clockwise around w, making one full rotation. Now consider the first 
component that is hit by r. As r rotates, the component hit changes exactly 
where there is a T- window of v. The changing sequence of components hit by r 
through one full rotation is the component sequence of v. This sequence starts 
and ends with Pg. 

It remains to show that the component sequence of i; is a Davenport-Schinzel 
sequence of order 2; since it is a sequence over h+1 symbols, this implies that the 
sequence length is at most 2h+l [26 , which in turn implies that there are at most 
2h T-windows of v. A Davenport-Schinzel sequence of order 2 is a sequence in 
which, for any two symbols a and b, the subsequence . . . , a, . . . , 6, . . . , a, . . . , &, . . . 
does not appear. 

The rest of the proof follows a simple geometric argument illustrated in Figure 
^ Assume that, in the rotation of r, it hits Ci at point ui, then hits Cj, then 
hits Ci at point U2 ■ Then the union of Ci and line segments vui and vu2, there 
is a Jordan curve that contains the relative interior of Cj in its interior region. 
Thus Cj will not be hit by r in its rotation after it passes U2- □ 



Lemma 3. For a polygon P with h holes, a fixed vertex v, and a fixed line 
segment £ in P, £ intersects at most 2(h + 1) windows of v: one right window, 
one left window, and 2h T-windows. 




Fig. 4. If, in rotational order, v 
sees Ci, then Cj, then d again, 
then it will never see Cj after 
that because Cj is 'trapped' in 
the interior of the red Jordan 



curve. 



Fig. 5. An illustration accompanying the proof 
of Lemma [3] The line segment £ is intersected by 
6 windows of v. Alternating invervals of visibility 
from V are shown in red and green. For £ to be 
intersected by more windows of v, more than 2 
holes would be required. 
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Proof. See Figure [sj Of the windows of p that cross £, at most one is a right 
window and at most one is a left window. To see this, consider a point p moving 
from one endpoint of £ to the other, directed so that Hne segment vp rotates 
clockwise around v as p moves, p can exit at most one left pocket of v and 
cannot enter a left pocket of u. p can enter at most one right pocket of v and 
cannot exit a right pocket of w. At most 2h T- windows oi p cross £ because, by 
Lemma |2] there are at most 2h T- windows oi p. □ 

Theorem 1. For a polygon P with h holes, 'Dy{P, V{P)) contains 0[{h + l)n'^) 
cells. 

Proof. By Lemma |3j each window is crossed by at most 2{h+l) windows of each 
vertex, and thus 2(ft, + l)n other windows in total. The total number of windows 
is ©(n^), so the number of points at which windows cross is 0{{h + The 
theorem follows by applying Euler's formula for planar graphs. □ 

2.3 Constructing and Using the Decomposition 

Chazelle and Edelsbrunner [7_ give an efficient algorithm for finding the pla- 
nar decomposition defined by a set C of line segments. The running time is 
C(|£|log|£| +k), where k is the number of intersections. In our case, |£| is 
C(n^) and Theorem [l] tells us that k is 0{(h + 1)71"^), so the total runtime for 
this step is 0{{h+l)n^). 

The visibility decomposition VyiP, V{P)) has an associated planar graph G. 
We are particularly interested in a customized graph dual G' of G that we can 
use to find cells of minimal visibility. To obtain G', we first build the dual graph 
of G. We then remove all edges corresponding to segments of dP. This gives us 
the underlying graph of G", and we will direct the edges as follows. 

Each edge of G that does not correspond to part of dP separates two cells; 
thus each edge of G' connects two cells. Consider an edge e of G' that connects 
two adjacent cells, or faces, fi and fj. If Vi and Vj are the respective sets of 
vertices that see points in fi and fj, then either Vi can be obtained from Vj by 
removing one vertex or vice versa. We direct e towards the face with the smaller 
corresponding set, i.e., towards the cell that is less visible. Doing this for all 
edges we obtain a directed acyclic graph. The sinks of this graph correspond to 
the cells of minimal visibility in 2?y(P, V{P)). 

Bose et al. used this technique and showed that, given the planar decompo- 
sition of a simple polygon P, in 0[n'^) time it is possible to construct G' and 
identify all sinks. This construction extends naturally to polygons with holes in 
time 0{{h+l)n^). 

It is worth noting that Bose et al. used the planar decomposition algorithm of 
Bentley and Ottmann [3] , which in this application is slower by a factor of log n. 
They do this because they claim the algorithm of Chazelle and Edelsbrunner [7] 
requires the line segments to be in general position (they will not necessarily be in 
general position even if the polygon's vertices are in general position). However, 
we have looked for and failed to find any original mention of this general position 
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Fig. 6. The visibility decomposition of a polygon with its directed dual edges indicated 
in red and its sinks, i.e., cells of minimal visibility, shaded green. 
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requirement. At any rate, degeneracies will only occur at the vertices of P, and 
an efficient workaround is possible. 

2.4 Bounding the Number of Sinks 



Lemma 4 ([5]). Given two right pockets pi and p2, of a point x, no point in 
Pi can see a point in p2. 

Proof. The proof of [Jl Lemma 8] works without modification when generalizing 
to polygons with holes. For completeness we restate it here. 

If a point a in pi could see a point b in P2, then the line segment between 
them must intersect both the window of pi and p2. Consider the window half 
planes of pi and p2- Since the two half planes intersect at x and both windows 
are right windows, one of the half planes must contain the windows of the other. 
The lemma follows from Lemma [TJ □ 

Lemma 5 ([5j). There is at most one point of intersection between all the right 
windows of a point x and all the right windows of a distinct point y. 

Proof. The proof of O Lemma 9] works without modification when generalizing 
to polygons with holes. For completeness we restate it here. 

Note that if a right window of x intersects a right window of y, then both x 
and y are visible from the intersection point. Also, from the fact that both are 
right windows, the base of one window must be contained in the pocket of the 
other. There are two cases to consider: either x is contained in a right pocket of 
2/ or a; is not contained in a right pocket of y. We start with the former. 

If x is contained in a right pocket of y, then by Lemma |4] x cannot see 
any other right window of y, and by Lemma [3] only one right window of x can 
intersect the right window of y's right pocket containing x. Therefore, the lemma 
follows in this case. 

Assume that x is not contained in a right pocket of y. Suppose that a right 
window ri of x intersects a right window r2 oi y. Since x is not conatained in 
a right pocket of y, the base of r2 must be contained in the pocket of ri. This 
implies that y is in the window half plane of ri . However, since y is visible from 
the intersection point of ri and r2, y must be in the pocket of ri by Lemma [TJ 
Therefore, the lemma follows since y is contained in a right pocket of x. □ 

Lemma 6. For distinct vertices Vi and Vj, if Vi is not in a right pocket of Vj 
then any T-window of Vi crosses at most one right window of Vj . 

Proof. By Lemma |4j a ray shot from Vi cannot leave a right pocket of Vj and 
enter another right pocket of vj. The lemma follows easily. □ 

Corollary 1. The total number of crossings between T-windows is 0{(h + l)^rt^) . 
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Theorem 2. For a polygon P with h holes, VyiP, V{P)) contains 0[{h + 
cells of minimal visibility, and this is tight in the worst case. 

Proof. Let Ri be the set of right windows of u^. We want to show that there are 
0{{h + l)n) sinks bordered by a window in Ri. We do this by showing that, for 
a vertex v, there are 0{h) sinks having an edge from a window in i?^, followed 
in clockwise order by an edge from a window of w. If u is in a right pocket of 
Vi, V can only see one window in Ri, namely the window bounding the pocket, 
and by Lemma |3] we are done. Otherwise, consider a sink s that has one edge 
formed by a right window r ^ Ri and the next edge in clockwise order formed by 
a window of v. The window of v must be either a right window or a T- window, 
otherwise s would not be a sink. By Lemma [Sj at most one right window of v 
can intersect a window in i?,. By Lemma [6j each T- window of v crosses at most 
one window in Ri. 

We have shown that there are 0{{h+ l)n) sinks bordered by a window in 
Ri. We can do the same for sinks bordered by windows in Li, the set of left 
windows of Vi. Thus the total number of sinks bordered by all right and left 
windows is 0((/i + l)n^). By Corollary [l] we know that the number of cells 
having two consecutive sides formed by T-windows is 0(^{h + l)^n'^), since each 
intersection point of two T-windows borders at most 4 cells. The upper bound 
of C((/i + for the total number of cells of minimal visibility follows. 

We prove the lower bound by example. See Figures [7j [8] [9] and[TOj 

□ 

3 Greedy Approximation 

Based on the previous section we treat the guarding problem as the abstracted 
problem of finding a minimum hitting set for a discrete range space S = {X, TV), 
where \X\ = n and \TZ\ = 0{{h + This range space can be obtained in 

time 0{{h+l)n^). 

Using the greedy algorithm for set cover/hitting set on this range space gives 
us a C'(logn)-approximation algorithm in OdXHT?-!) — 0{{h + Vfn^) time. 

4 Improved Approximation via e-Nets 

We now turn our attention to the task of achieving a better approximation 
factor when opt is small; for our bounds on time complexity we can assume 
that OPT = 0(n}/^^. Otherwise an approximation factor of 0(\ogn) is also 
©(log opt). 

4.1 VC-Dimension and e-Nets 

We use standard techniques for approximating hitting set that are based on the 
concept of VC-dimension first introduced by Vapnik and Chervonenkis |28j in 
the area of learning theory. 
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Fig. 7. A polygon with 0{{h + l)^n^) sinks. Eacii light strip, indicated in yellow, is the 
intersection of visibility polygons of two adjacent vertices. This polygon has 10 = &(n) 
light strips and 6 holes. The construction easily generalizes to higher values of n and 
h. This figure is diagrammatic; in the real polygon the holes would be closer to the 
region in which horizontal and vertical light strips interact. 




Fig. 8. Detail of one of the light strips. The shadow strips, indicated in light grey, are 
the regions not seen by the vertex at the bottom right. Each light strip has h/2 shadow 
strips. 
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Fig. 9. Detail of tlie interaction of 
shadow strips. Tiiere are ©((/i + l)^n^'j 
dark grey squares, eacli of wliicli is the 
intersection of shadow strips and must 
contain a sink. 



Fig. 10. Detail of the intersection of 
shadow strips. Dual edges are indicated 
in red; since they are all pointing in- 
wards, the dark grey region must con- 
tain a sink. 



Definition 1 (VC-Dimension |28p . For a range space S = (X.TZ), let Y be 
a maximum cardinality subset of X such that TZr\Y^2^. The VC-dimension 
of S is equal to \Y\ . 

We use the following bounds on the VC-dimension of visibility systems in 
polygons due to Valtr '27\. 

Theorem 3 (|27j). The visibility system of a polygon with h holes has VC- 
dimension at most 23 if h — (i.e., if the polygon is simple), and 2 log2 h + 
41og2log2/i + o(l) =0{l + logh) ifh>l. 

In systems of bounded VC-dimension, small e-nets can be constructed via 
random sampling. We use the following result of Blumer et al. f?I . 

Theorem 4 (|4])- Por a measure fi on the elements of a range .space of VC- 
dimension d, and for any e G (0,1] and any 5 G (0,1], a random sample of 
m{e, 5) elements drawn according to fi forms an e-net with probability at least 
1 — 6, where 

(4 2 8d 13\ 

m(e, 0) — max - log - , — log — 
\e e e J 

The following is a straightforward consequence of Theorems [3] and |4] 

Corollary 2. There exists a function f{h,e) — O^^l -f log(/i -I- 1)^ Mog 

such that, for a measure pL on a polygon with h > holes, and for any e G (0, 1], 
a random sample of f{h, e) points drawn according to /i forms an e-net with 
probability at least 1/2. 
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4.2 Approximation via e-Nets 

Bronnimann and Goodrich [5] developed a method for turning algorithms for 
finding £-nets into approximation algorithms for finding minimum hitting sets. 
The key is in finding an optimum, or approximately optimum, measure /x on 
the elements of a range space. Their algorithm essentially 'learns' a measure 
/Li through an iterative doubling technique. The optimum measure ^* is the 
distribution that maximizes the value of e* such that every e*-net is a hitting 
set. In fact, the size of a minimum fractional hitting set is exactly OPT/ — 1/e* 
[12]; this serves as a lower bound for opt. 

Iterative doubling. The B&G algorithm [B] §3.1] finds a measure such that 
every e'-net is a hitting set for some e' > ^ Jprj. . The algorithm starts by assign- 
ing a weight of 1 to every element (when the algorithm terminates these weights 
are normalized to obtain n'). In each iteration the algorithm makes one call to 
the e-net finder, which returns an £-net Y Q X and one call to a verifier. The 
verifier checks if the given set F is a hitting set. If so, the algorithm returns Y 
and terminates. If not, the verifier returns a range R that is not hit by Y; the 
algorithm doubles the weight of every element in f?, then starts a new iteration. 
The algorithm is guaranteed to terminate after ^^(opt • log \X\) iterations and 
the total weight of elements cannot exceed \X\'^ [B]. Since the B&G algorithm 
does not know the value of opt or e' a priori, it must make several guesses, 
starting at a constant value such as 1/2 and halving the guess after each failed 
run of the algorithm. The result of this is that the entire algorithm must be run 
©(log opt) times 0. 

Random sampling and verification. We show that, for our application, the total 
time complexity of random sampling and verification is ©(n'^). 

For verification we build a directed bipartite graph with n vertices represent- 
ing elements and 0(n^^ elements representing ranges, where an element vertex 
is adjacent to a range vertex if and only if the element is in the range. This 
graph is constructed in ©(n"^) time and can be constructed once and used for all 
iterations of the algorithm. Each range vertex is given a boolean fiag indicating 
whether or not it is hit by the input set. In each verification round we do the 
following: 

1. In 0(1 7?. I) time, reset the flags. 

2. In C(|l^| • \Tl\) total time, for each element in the input set Y, mark all 
incident ranges as 'hit'. 

3. In 0(|7?,|) time, scan through the ranges to find an imhit range or verify that 
all ranges are hit. 

4. In 0(|7?,|) time, if there is an unhit range, double the weight of every element 
in that range. 

Thus each verification round takes C(|i^| • \Tl\) time. 

For random sampling, we assume that we can sample a random bit in 0(1) 
time, allowing us to sample uniformly from an array of k elements in C(log k) 
expected time. Through all iteratoins of the algorithm, every element weight is a 
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power of 2. Since the total element weight does not exceed \X\'^, each element has 
its weight doubled C(log \X\) times in each run of the B&G algorithm. We can 
maintain a partition of the elements based on how many times each element's 
weight has been doubled. In this way we can sample an element according to 
the weight function in ©(polylog |X|) amortized time. The additional cost of 
maintaining this partition is ©(IXj log \X\) time per run of the algorithm. 

Total time complexity. For our application, we have that Y = 0(n^^^ polylog n) 
and \TZ\ = 0{{h+ In each iteration of the B&G algorithm, the time 

complexity of sampling for the e-net finder is 

f{h, e)0(polylog n)=o((l + \og{h + 1)) 1 log - polylog n 



e 

— 0(0PT • log OPT • polylog n) 
= ©(^n^/^ polylog . 

The time complexity of the verifier is 0(1^1 • \7l\) — 0[{h + l)^n''/'^ polylogn). 
We perform ©(log opt) runs of the algorithm and in each run there are 
C(OPT • logn) iterations. Thus the total number of iterations is ©(opt polylogn) 
= 0(n^^^ polylog nj . Our total running time, including all runs of the algorithm, 
all iterations of the e-net finder and verifier, and additional overhead, does not 
exceed 0{{h + l)^n^). 

Approximation ratio. When the last run of the algorithm terminates with suc- 
cess, we are left with a measure /i' on X, along with an e'-net that is a hitting 
set, where e' > ^Tfjpji- The size of this hitting set is at most 



o(^(l + log(/i+l))llog^ 

0( ( 1 + log(/l + 1) JOPT • log OPT 



Thus we have shown the following theorem. 

Theorem 5. For simple polygons or polygons with the number of holes bounded 
by a constant, there exists an approximation algorithm running in 0(n^^ time 
with an approximation ratio o/ ©(log opt). For polygons with h>2 holes, there 
exists an approximation algorithm running in 0{h^n^^ time with an approxima- 
tion ratio o/ ©(log log opt). 



5 Further Improved Approximation for Simple Polygons 

For simple polygons. King and Kirkpatrick |19) presented an e-net finder that 
returns e-nets of size ©(Moglog j); they did not analyze its time complexity, 
rather they simply stated that it runs in polynomial time. Here we show that it 
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runs in time 0(n^ log log ^). This means that, for simple polygons, we can use 
it to replace the random sampling e-nct finder of the previous section without 
pushing the total running time above 0(n'^). 

First we note that the B&G algorithm always uses values of e that are powers 
of 2. This means that 

which slightly facilitates analysis. With this restriction on e, the K&K net finder 

recursively partitions the vertex set of the polygon, keeping its cyclic ordering. 
Define t = [log log t\ . At the i^^ level of the partition, each subset is further 
divided into hi subsets, with 

_ /22*-'+i.4f-2i-* ,i = l 
'~\22*-*+i ,l<i<t. 

If fi is the number of new fragments created by the i*^ fragmentation step, this 
gives us 

r 1 ,i = o 

fi=lAt- 22*-2'--t+i+i , < i < t 
I 4t • 22* ,i = t. 

In O(nloglogi) time we can build a tree corresponding to the hierarchical 
decomposition in which each node stores the corresponding subset in cyclic order. 
This tree has loglog j) nodes. Determining guards to place for a pair of 
sibling fragments, respectively storing vertex sets Ui and U2, can be done with 
0(|U'i||t/2|) calls to a visibility oracle. The visibility matrix of the vertices of P 
can easily be built in 0{n^) time and serves as a constant-time oracle. 

Each of the nodes at level z — 1 contains 0{n/ fi-i) vertices; it has hi 
normal children, each with 0{n/fi) vertices, and one dummy child with 0{n) 
vertices. The cost of placing all guards at level i is therefore 




Thus the cost of placing all guards at all levels is O^ti^ log log j), as desired. 

Using this e-net finder instead of random sampling, we can achieve an ap- 
proximation ratio of O(loglogOPT) for simple polygons in O(n^) time. 

Theorem 6. For simple polygons, there exists an approximation algorithm run- 
ning in O(n^) time with an approximation ratio 0/ C?(log log opt) . 
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